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ABSTRACT 

We report the discovery in space of a disilicon species, SiCSi, from observations 
between 80 and 350 GHz with the IRA1V0 30m radio telescope. Owing to the close 
coordination between laboratory experiments and astrophysics, 112 lines have now 
been detected in the carbon-rich star CW Leo. The derived frequencies yield improved 
rotational and centrifugal distortion constants up to sixth order. From the line profiles 
and interferometric maps with the Submillimeter Arra>12[ the bulk of the SiCSi emis¬ 
sion arises from a region of 6" in radius. The derived abundance is comparable to that 
of SiC 2 . As expected from chemical equilibrium calculations, SiCSi and SiC 2 are the 
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most abundant species harboring a Si-C bond in the dust formation zone and certainly 
both play a key role in the formation of SiC dust grains. 

Subject headings: Stars: individual (IRC +10216) — stars: carbon — astrochemistry 
— stars: AGB and post-AGB 


1. Introduction 


Interstellar dust grains are synthesized in two main types of sources: the inner winds of 
asymptotic giant branch (AGB) stars and the ejecta of massive stars, mainly supernovae (SNe). 
Dust grains are formed from molecular seeds. In AGB stars, molecules such as TiO, VO, ZrO, 
C2, CN, and C3, among others, are known to be present in their photospheres since the beginning 
of the 20th century. Because of the high stability of the CO molecule, depending on whether the 
C/O ratio is >1 (C stars) or <1 (M stars), the gas becomes either carbon- or oxygen-rich after CO 
reaches its equilibrium abundance, leading to a different chemistry and therefore to different dust 
formation schemes. In SNe, different types of dust can be formed depending on the degree of 
enrichment in heavy elements. The dust formation efficiency in SNe is still a matter of debate. Re¬ 
cent studies with the Herschel Space Telescope point to much larger dust masses than previously 
thought (iMatsuura et al. 201 lj, 2015 ). 


Dust formation can be simplified as a two-step process: formation of nucleation seeds fol¬ 
lowed by grain growth through condensation of refractory molecules at high temperature and other 
less refractory species at larger distances from the star. There are, however, many mysteries in 
this picture of events, starting from the fundamental step of the formation of the nucleation seeds, 
which essentially are refractory species (Gail 120100 . The presence of SiC grains in C-rich AGBs 
was confirmed by the detection of an emission band at ~11.3 /urn (Treffers & Cohen 1974). This 
band has been found towards a large number of C-rich stars with the IRAS and ISO satellites. 
However, the molecular precursors of SiC dust grains ar e still unknown. SiC molecules have bee n 
detected in the external shells (>300 R*) of IRC +10216 ( Cemicharo et al. 1989 : Patel et al. 2013 ). 


In the inner layers, the most abundant molecule harbouring a Si-C bond is SiC 2 (IThaddeus et al. 


1984 ICernicharo et al.l l2010h . whereas the most abundant Si-bearing species is SiS, which, to¬ 


gether with SiO and SiC 2 , account for a significant fraction of the available silicon (lAgundez et al. 

20121 ) . 


CW Leo, located ^ 130 pc from us^is a Mira variable star with a period of 630-670 days and an 
amplitude of -1 mag in the K band (Menten et al. 2012). It is one of the brightest infrared sources 
in the sky. Due to its close proximity, IRC +10216, the circumstellar envelope of CW Leo, has 
attracted many studies because it is exceptionally rich in molecular species. Half of the known 
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interstellar species are observed in this C-rich envelope. The observed molecules range from 


Cemicharo et al. 2000, 2010), to molecules containing refractory elements ( 

Cemicharo & Guelin 


1987) 

, and long carbon chain species C,,H and their anions ( 

Cernic 

laro & Guelin 1996; 

Cemicharo et al. 

OO 

O 

O 

<N 

Guelin et al. 1997; McCarthy et al. 2006; Thaddeus et al. 

2008 

, and references therein). 


Among the species detected in this source, the silicon-carbon species Si„C,„ could play an impor¬ 


tant r ole as seeds of SiC dust grain s. The simplest members of such a family, SiC (ICernicharo et al. 


1989b . and SiC 2 (Th ad deus et al. 19841). are known to be present in IRC +10216. However, Si 2 C 


hereafter SiCSi), which is predicted to be very abu ndant from chemical equilibrium calculations 
Teiero & Cemicharo 1991; Yasuda & Kozasa 2012) and could play a key role in the formation of 
SiC dust grains, has not been found so far. 


In this Letter we present the detection of SiCSi in IRC +10216 based on a close coordination 
between laboratory, astrophysics and radio astronomy. This synergy has resulted in the detection 
of 112 lines of this new species using data from the IRAM 30-m radio telescope. Nine of these 
lines are also detected with the Submillimeter Array (SMA). 


2. Observations 


In the course of searches for new molecules we have covered a large fraction of the 3, 2, 1 
and 0.8 mm spectrum of IRC+10216 with a high sensitivity using the 30-m IRAM radio telescope. 
In the 3 mm window, the data acquired during the last 30 years cover the 80- IT 6 GHz domain 
with very high sensitivity (1-3 mK) . Examples of thes e data can b e found i mAgundez et al.1 (120081 


2014) . ICernicharo & Guelinl (11996b . ICernicharo et al.1 (120071. l2008h . and references therein. Most 


of the 2 mm data come from the line survey of IRC +10216 carried out by Cemicharo et al. ( 20001) . 
com plemented with additional data obtained during the search for spe cific molecular species (see, 
e.g. Guelin et aklhoQdl: lAgundez et al.ll2008[ I 2 OI 2 I: iFonfrfa et al.l 20061) . The sensitivity of the 2 mm 
observations varies between 1.5 and 10 mK. The 1 mm and 0.8 mm data come from observations 
carried out during the searches quoted above and a line survey between 290 and 355 GHz carried 
out using the new EMIR receivers in 2010 and 2011 with a sensitivity of 2-5 mK, depending on 
the atmospheric transmission. 


For frequencies above 250 GHz the spectrometers were two autocorrelators with 2 MHz of 
spectral resolution and 4 GHz of bandwidth. For all other observations the spectral resolution was 
1 MHz provided by filter banks or autocorrelators. The observing mode, in which we wobbled the 
secondary mirror by +90” at a rate of 0.5 Hz, and the dry weather conditions (sky opacity at 225 
GHz was below 0.1 in most observations) ensured flat baselines and low system noise temperatures 
(T sys - 100-400 K depending on the frequency). 
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Fig. 1.— Selected transitions of SiCSi among the 112 observed lines. Unidentified (U) lines are 
indicated by their frequency in MHz. The intensity scale is antenna temperature in mK. The data 
are shown as histograms. Violet continuum lines show the predicted line profile from the best 
model for the abundance of SiCSi (see text). 
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The selected observing method, with the off position located at 180" from the star, provides 
reference data free from emission for all molecular species but CO (see Cernicharo et al.l [20150. 


The emission of all other molecular species is restricted to a region < 20-30" (see, e.g., Guelin et al. 


19931) . The intensity sc ale, antenna temperature (T*), was corrected for atmospheric absorption 
using the ATM package ( Cemicharolll985 : Pardo et aPhoQl ). The main beam antenna temperature 
can be obtained by dividing T* A by the main beam efficiency of the telescope which is 0.81, 0.59, 
and 0.35 at 86, 230, and 340 GHz, respectively. Calibration uncertainties for data covering such 
a large observing period have been adopted to be 10%, 15%, 20%, and 30% at 3, 2, 1, and 0.8 
mm, respectively. Additional uncertainties could arise from the line intensity fluctuation with 
time induced by the variation of the stellar infrared flux, which has been recently discovered by 
Cernicharo et al. (2014). All data have been analyzed using the GILDAS package^]. 


The data revealed several hundreds of spectral lines which could not be assigned to any known 
molecular species collected in the public CDMS (M u lle r et al. 2 005) and JPL (Pickett 1998) spec¬ 
tral databases and in the MADEX code ( Cernicharo 12012 ). Most of these lines show the character¬ 
istic U-shaped or flatted profiles with linewidths of 29 km s _1 . However, above 250 GHz a signifi- 
cant number of l i nes are very narrow and come from the dust formation zone of CW Leo (see, e.g., 
Patel et al.l 1201 ll: ICernicharo et al.l 120130 . Among the unidentified lines observed with the IRAM 
30-m telescope, we have been able to assign 112 to the rotational spectrum of SiCSi through an 
iterative procedure described below and a close synergy between molecular spectroscopy and as¬ 
trophysics. Some selected lines among those observed in IRC +10216 are shown in Figure 1. Nine 
lines observed above 300 GHz with the IRAM 30-m telesc opes a nd assigned to SiCSi correspond 
to unidentified lines observed with the SMA by Patel et al. (20IT). The spatial distribution of some 
of these lines, as derived with the SMA, is shown in Figure 2. 


3. Spectroscopic constants of SiCSi 


Disilicon carbide (SiCSi) has a Ci_ v symmetry and a '/ 

ii electronic ground state with a modest 

permanent dipole moment of ~1 D along the h-inertial axis 

(Gabriel et al. 1992; 

Barone et al. 1992; 

Bolton et al. 

1992; 

Spieffiedel et al. 

1996; 

McCarthy et al. 

2015 

). Because the two equivalent off- 


axis silicon atoms are bosons, only half of the rotational levels exist ( K a + K, even). 


The search for hnes of SiCSi in space began with the lines measured in the laboratory by 


M cC arth y et ah (2015b . which allowed us to accurately predict frequencies with K u - 0,1 in the 


3 mm domain. Although the molecule is predicted to be rather floppy, with deviations up to 10- 


1 http://w w w. iram. fr/IR AMFR/ GILDAS 
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Fig. 2.— Maps of the average of four SiCSi lines: 183.15 - I82.1A, I63.11 - 162,14, 14 3 ,n - 14 2 .i2 
and 12 3 9 - 12 2 ,io (from the SMA line survey of IRC +10216 Patel et al. (2011)1. Their visibilities 
were combined before imaging. The upper left panel is the integrated intensity map, while the 
remaining ones are the channel maps with Vlsr velocities indicated on top right in km s _1 . The 
size of the synthesized beam is shown in bottom right panel. The halftone image in the first panel 
is scaled linearly from 0 to 2.7 Jy beam” 1 km s _1 and the contour levels (in the same units) have 
the starting value and intervals of 0.4 in the first panel and 0.04 in remaining panels. 
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20 MHz for transitions higher in J than those measured in the laboratory (J < 10), we quickly 
found 3 lines in the 3 mm spectrum of IRC +10216 that could be analyzed simultaneously with the 
laboratory lines. From the newly derived spectroscopic constants, which now include additional 
centrifugal distortion terms, we searched for additional lines until we found all the intense lines 
predicted in the 3 mm domain with upper level energies below 100 K. Using the spectroscopic 
constants derived from a combined fit to the laboratory and IRC +10216 data, we continued the 
search for K u - 0,1 lines in the 2 mm and 1 mm domains, and up to 20 additional lines were 
easily identified. Differences between the predicted and observed frequencies were again 10-50 
MHz, but because the lines are so intense (typically 10-20 mK) and have a distinctive line shape, it 
was not difficult to make additional assignments (see Figure 1). Nevertheless, on occasion, we had 
to discard several initial assignments because subsequent predictions failed to predict new lines. 
With the laboratory lines and the ~30 astronomical lines involving K a = 0 and 1 levels, a refined 
set of constants were derived, and we then proceeded to search for lines with K a = 2,3 in the entire 
30 m data set of IRC +10216 that was available to us. With these new astronomical measurements 
in hand, several additional lines with frequencies up to 180 GHz were measured in the laboratory. 
Most importantly, several of these lines correspond to transitions involving K a = 1 and 2 levels, 
which enabled us to precisely determine the D K term, and confirm the astronomical assignments 
with higher K a . 


Ultimately, we were able to assign 112 lines in IRC +10216 to SiCSi with / < 48 and K a < 5. 
At least 30% are free of blending, and the uncertainties of the derived frequencies are less than 
1 MHz. Other lines were slightly blended, but from their characteristic line profile it was still 
possible to measure their frequencies with an accuracy of ~1 MHz. Above 250 GHz, most ob¬ 
servations have a spectral resolution of 2 MHz and the frequencies are accurate to 1-1.5 MHz. 
From the unblended lines, an average expansion velocity of 14.0+1.0 km s' 1 was derived, which 
is similar to that obtained for most lines in IRC +10216 (ICernicharo et al.1 20001) . Hence, the dis¬ 
covery of SiCSi in IRC +10216 -which is based on 112 observed lines- is one of the most robust 
identifications ever published for a first detection of a molecule in space. Some 20 additional lines 
have been also identified, but they are so heavily blended that accurate frequencies and intensities 
could not be derived. 


The rotational constants were derived with a Watson Hamiltonian that included distortion con¬ 
stants up to sixth order in the A and S reductions (Table 1). The fit with the S reduction was done 
with the SPFIT program Pickett ( 19911) . while that in the A reduction was done with fitting pro¬ 
grams associated with MADEX (Cerni charo 2012). We confirmed that both SPFIT and MADEX 
yielded the same results in the A reduction. The rotational constants in the A and S reductions 
agree rather well, but the S reduction is preferred because S K and h JK are three orders of magni¬ 
tude larger than the corresponding constants d\ and h 2 and the predicted frequencies for high J or 
K a will be less accurate. Because there are only a few observed lines with K a > 3, further improve- 


















Table 1: Spectroscopic constants (MHz) of SiCSi in comparison to SiC 2 


Constant 

S-reduction 

A-Reduction 

SiC 2 A-reduction 

A 

64074.3366(44) 

64074.33623(37) 

52474.1930(677) 

B 

4395.51772(41) 

4395.621072(844) 

13158.71537(204) 

C 

4102.13098(62) 

4102.02789(107) 

10441.58029(168) 

A j xlO 3 

9.66776(224) 

9.73150(178) 

13.20278(397) 

A j K 

-0.856833(73) 

-0.8572075(610) 

1.5382238(882) 

A K 

23.58788(178) 

23.58805(148) 

-1.2159(174) 

d\/6j x 10 3 

-1.52630(34) 

1.519832(437) 

2.41191(253) 

di/S k xl(h 

-0.00318(32) 

5.1591(454) 

8.70880(415) 

Hj x 10 s 

-3.627(116) 

-4.1349(949) 

-8.312(460) 

Hjk x 10 5 

1.9373(67) 

1.93298(547) 

-5.0158(993) 

H kj x 10 3 

-1.8882(95) 

-1.88755(791) 

0.39025(399) 

H k x 10 2 

4.8632(204) 

4.8633(168) 

0.2487(529) 

h\/hj x 10 9 

-5.304(225) 

-5.231(187) 

-3.46(271) 

h 2 /h JK x 10 9 
h 3 /h K x 10 3 

-2.614(255) 

-6586(361) 

-35151(191) 
1.1162 (286) 


Note. — 1 cr uncertainties (in parentheses) are in the units of the last significant digits. 
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ment of the rotational constants was not possible. We estimate that the calculated frequencies with 
J < 60, K a < 5 are reliable up to 400 GHz. The list of observed lines in IRC +10216, together 
with the derived line parameters, energies of the upper levels, and line strengths is provided in 
Table 2 (online), which also gives the frequencies and uncertainties for the 22 lines observed in the 
laboratory ([McCarthy et al. 2015 ). 


The SiCSi distortion constants are as large as those of SiC 2 , a well-studied molecule also 
considered to be very floppy. Table 1 provides, for comparison purposes, the rotational constants 
derived from a fit to the lines of SiC 2 reported by Muller et al. (2012) including the same constants 
as for SiCSi plus additional higher distortion constants. Both, A K and H K are one order of mag¬ 
nitude larger for SiCSi than for SiC 2 , while A JK and d K are of the same order despite SiCSi being 
somewhat heavier than SiC 2 . Both molecules have a very low-lying vibrational mode, around 100- 
200 cm" 1 . Hence, by analogy with SiC 2 for which the antisymmetric v 3 mode is readily observed 
in IRC +10216, we might also expect to observe rotational lines of SiCSi in its symmetric mode 
v 2 . Nevertheless, SiCSi has its dipole moment along the /?-axis while in SiC 2 it is aligned along the 
a-axis which produce a significantly different rotational spectrum. 


4. Discussion 


The observed line profiles of SiCSi go from U-shaped, in the case of low excitation lines lying 
at high frequencies and thus observed with a smaller beam, to flat-topped, for high-excitation lines 
lying at low frequencies (see Figure 1). Such behaviour indicates that SiCSi is concentrated around 
the star but with a relatively extended distribut ion, as occurs in the case of metal-bearing species 
(ICemicharo & Cue1inlll987l:ICuelin et al .III 9931) . The SMA maps, which correspond to transitions 
involving upper level energies between ~60 and ~100K, show that the emission is concentrated 
in a region of - 6" in diameter (see Figure 2). A slightly larger brightness distribution size can 
be expected for lines involving lower energy levels (see Figure 1). We have checked if SiCSi 


could be responsible for some of the lines detected with ALMA in IRC +10216 (Cernicharo et al. 


2013]). Unfortunately, only three strong lines are within the spectral domain covered by these 


ALMA cycleO observations and only one seems to be free of blending and is clearly detected in 
the central pixel. However, these observations were selected to trace the dust formation zone and 
their sensitivity is not good enough to trace the spatial distribution of weak lines. 


To initially estimate the abundance and origin of the SiCSi emission from transitions involving 
different energies, we show in the top panel of Figure 3 a rotational diagram based on the 112 
observed lines. We have adopted a gaussian source size with a radius 5". The observed integrated 
intensities have been corrected for dilution in the beam and for the main beam efficiency of the 
telescope. We can distinguish three zones with different rotational temperatures. The lines with 




















Fig. 3.— The top panel shows the rotational diagram of the 112 observed lines of SiCSi. The 
three zones discussed in the text are indicated by vertical dashed lines. The bottom panel shows 
the radial abundance distribution of SiCSi, SiC 2 , and SiC as calculated by chemical equilibrium 
(dashed lines) and as derived from the observations of IRC+10216 (solid lines). The derived 
abundance profile of SiC 2 is taken from ICemicharo et al. (2010). 
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upper level energies below 60 K (zone I) have been fit with a rotational temperature, T rot , of 
31 + 5 K and a column density of (1.4 + 0.4) x 10 15 cm' 2 . The lines with upper level energies 
between 60 and 200 K (zone II) have been fit with T rot = 118 + 12 K and a column density of 
(3 ± 0.7) x 10 15 cm -2 . Finally, the lines with upper level energies above 200 K (zone III) have 


been fit with T rot = 280 + 50 K and a column density of (3.2 ± 0.9) x 10 13 cm . The rotational 


temperature derived for zone 
r ^14" in SiC 2 (Guelin et ah 


could correspond well with the extended molecular ring observed at 
1993). Its value agrees well with the rotational temperature derived 


for other species (lAgundez et al.l 120081: ICemicharo & Gueli nil 19961) . Zone II corresponds to the 
emission region of metal-bearing species and high energy lines of SiC 2 with r < 5" (Figure 2 and 
Velilla Prieto et al.lhpl 5 ). Finally, the high energy lines (zone III) correspond to the region with r < 
2" trace d by high-7 lines and yibrationally excited states of Si C 2 , SiS, HCN, H NC among other 


specie s (ICemicharo et al.ll2010L 1201 ll: I Velilla Prieto et al.ll2015l) . In this region, ICemicharo et al. 


( 2010() derived for SiC 2 T rot -204 K and a column density of ~ 8 x 10 cm* 2 averaged over a 
beam of 30". When corrected for the source radius of 5" adopted in this work, the column density 
of SiC 2 becomes 3 x 10 15 cm* 2 . Hence, SiCSi is as abundant as SiC 2 in the dust formation zone 
of IRC +10216. The large difference in the intensities of the observed lines of SiCSi and SiC 2 , 
near a factor of 100, arise from the lower dipole moment and the larger partition function of SiCSi. 
When corrected for these effects, the observed intensities of the individual lines of SiCSi and SiC 2 
indicate a SiC 2 /SiCSi abundance ratio of about 5. 


Early chemical equilibrium calculations by Tejero & Cernicharo (1991) showed that in the 
inner envelope of C-rich AGB stars the most abundant gas phase species containing a Si-C bond 
are SiC 2 and SiCSi^while SiC is predicted to be present with a much lower abundance (see also 
Yasu da & Kozasa 2012 and Figure 3). To constrain the abundance and spatial distribution of SiCSi 
in IRC +10216 we have carried out radiative transfer calculations to obtain the emergent line pro¬ 
files^ and compare them with the observed ones. The physical model of the envelope is taken from 
Agundez et al.l (120121) with the downward revision of the density of particles in the regions inner 
to 5 R* derived by I Cem ic har o et al.l (120131) . Since collisional rate coefficients are not available for 
SiCSi we have assumed local thermodynamic equilibrium, which is a reasonable approximation 
given the low dipole moment of SiCSi. We find that the intensities and profiles of the low ex¬ 
citation lines observed can be well reproduced by adopting a SiCSi abundance relative to H 2 of 
4 x 10“ 8 , from the inner layers out to the photodissociation region. We have assumed that SiCSi is 
photodissociated by the ambient interstellar ultraviolet field with a rate of 10“ 9 x exp (-2 A v ) s* 1 , 
where A v is the visual extinction. However, to reproduce the high intensities observed for lines 
involving upper level energies above 200 K, it is necessary to increase the abundance of SiCSi in 
the warm inner regions to 2 x 10* 7 (see Figure 3). Finally, a depletion in the very inner regions, 
which is consistent with the expectations from chemical equilibrium, yields a better agreement 
with the observed line profiles (see Figure 1). The derived abundance profile for SiCSi is shown in 
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Figure 3. 

High angular resolution interferometric observations are needed to better constrain the abun¬ 
dance gradient of SiCSi in the inner layers. In any case, we find that disilicon carbide is as abundant 
as silicon dicarbide in the inner layers. The depletion in the abundance of SiCSi at 50 R*, which 
leads to a SiC 2 /SiCSi abundance ratio of ^5 in the outer regions, would be consistent with its 
incorporation into SiC dust grains. 

We thank Spanish MINECO for funding under grants AYA2009-07304, AYA2012-32032, 
CSD2009-00038, and ERC under ERC-2013-SyG, G.A. 610256 NANOCOSMOS. The new labo¬ 
ratory measurements in Cambridge were supported by NASA Grant NNX13AE59G. 
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Table 2: Observed lines of SiCSi in the laboratory and in IRC +10216 


a 

Ka 

Kc)„ 

<,j 

K a 

Kc)i 

v^ ? (unc) 

MHz 

jTjxrfv 

K kms 1 

Predicted^+^o (unc) 
MHz 

K 

S 

cr 

mK 

5 

1 

5 

6 

0 

6 

6718.9927 (0.002) 


6718.994 (0.001) 

8.9 

2.586 

lab 

8 

0 

8 

7 

1 

7 

12018.3548 (0.002) 


12018.353 (0.001) 

14.7 

3.692 

lab 

20 

1 

19 

19 

2 

18 

19095.4772 (0.002) 


19095.477 (0.002) 

89.8 

4.991 

lab 

15 

2 

14 

16 

1 

15 

24506.4805 (0.002) 


24506.481 (0.002) 

60.4 

3.644 

lab 

3 

1 

3 

4 

0 

4 

24959.6919 (0.002) 


24959.692 (0.001) 

5.3 

1.528 

lab 

10 

0 

10 

9 

1 

9 

31198.7591 (0.002) 


31198.761 (0.002) 

22.4 

4.861 

lab 

1 

1 

1 

2 

0 

2 

42663.0910 (0.002) 


42663.088 (0.002) 

3.3 

0.504 

lab 

2 

1 

1 

2 

0 

2 

60248.2200 (0.020) 


60248.209 (0.003) 

4.1 

2.493 

lab 

4 

1 

3 

4 

0 

4 

61298.7200 (0.020) 


61298.746 (0.003) 

7.0 

4.450 

lab 

6 

1 

5 

6 

0 

6 

62974.9200 (0.020) 


62974.923 (0.005) 

11.6 

6.340 

lab 

8 

1 

7 

8 

0 

8 

65310.0300 (0.020) 


65310.040 (0.008) 

17.8 

8.135 

lab 

1 

1 

1 

0 

0 

0 

68154.6000 (0.020) 


68154.606 (0.004) 

3.3 

1.000 

lab 

10 

1 

9 

10 

0 

10 

68348.4600 (0.020) 


68348.474 (0.011) 

25.7 

9.806 

lab 

3 

1 

3 

2 

0 

2 

84424.5800 (0.020) 


84424.602 (0.006) 

5.3 

2.002 

lab 

5 

1 

5 

4 

0 

4 

100120.6600 (0.040) 


100120.678 (0.009) 

8.9 

3.022 

lab 

8 

2 

6 

9 

1 

9 

109709.6000 (0.040) 


109709.645 (0.018) 

26.2 

1.447 

lab 

7 

1 

7 

6 

0 

6 

115262.8500 (0.040) 


115262.899 (0.013) 

14.1 

4.073 

lab 

11 

1 

11 

10 

0 

10 

144033.4200 (0.040) 


144033.475 (0.023) 

29.3 

6.327 

lab 

3 

2 

2 

3 

1 

3 

180036.9000 (0.040) 


180036.916 (0.017) 

13.9 

1.456 

lab 

5 

2 

4 

5 

1 

5 

181405.4500 (0.040) 


181405.432 (0.014) 

17.6 

2.548 

lab 

7 

2 

6 

7 

1 

7 

183384.5800 (0.040) 


183384.603 (0.012) 

22.9 

3.561 

lab 

9 

2 

8 

9 

1 

9 

185976.7800 (0.040) 


185976.772 (0.012) 

29.8 

4.527 

lab 

16 

1 

15 

16 

0 

16 

82260.5000 (0.500) 

0.109 (22) 

82260.891 (0.029) 

59.2 

13.826 

1.7 

3 

1 

3 

2 

0 

2 

84424.1900 (0.330) 

0.107(18) 

84424.602 (0.005) 

5.3 

2.002 

1.4 

9 

2 

8 

10 

1 

9 

86429.8400 (0.500) 

0.059 (10) 

86429.537 (0.008) 

29.8 

1.888 

1.1 

18 

1 

17 

18 

0 

18 

88716.0400(1.150) 

0.091 (14) 

88715.413(0.040) 

73.7 

14.774 

1.7 


Note. — 1 cr uncertainties (in parentheses) are in the units of the last significant digits. The full table containing all 
observed lines can be found in the online version of the paper. 
















